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The behavior of polyelectrolytes and polyampholytes in semi-dilute solutions is investigated 
theoretically. Various statistical charge distributions along the polyelectrolyte chains are considered: 
(i) smeared, where the charges are uniformly distributed along the chain, (ii) Annealed, where the 
charges are allowed to associate and dissociate from the chain, (iii) Permuted, where the total 
number of charges on the chain is fixed, but the charges can move along the chain, (iv) Quenched, 
oo ■ where the charges on the chains are "frozen" in a random configuration. Finally, we also consider 

(v) polyampholytes, where each monomer can be positively or negatively charged, or neutral. Path 
integral formulation was used to derive mean field free energies for the different models. Self- 
consistent field equation is obtained for the polymer order parameter and a Poisson-Boltzmann like 
equation for the electrostatic potential. We show that the difference between the permuted and the 
C/3 ' smeared models is a constant shift in the chemical potential leading to similar mean field equations. 

Within mean-field the quenched model is found to be equivalent to the annealed one, provided that 
the system is coupled to a reservoir of polyelectrolyte chains. The random phase approximation is 
used to calculate the monomer-monomer structure factor S(q) for the different statistical charge 
distribution models. We show that in the annealed model, fluctuations of the the monomer charges 
contribute to the electrostatic screening in addition to the free ions in the solution. The strength 
of this screening depends on the variance of the monomer charge distribution and is especially 
important for polyampholytes in bad solvent conditions where the mesophase separation is enhanced. 
The ratio between the variance and the net average charge determines whether polyampholytes 
behave as polyelectrolytes or as neutral chains. 
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I. INTRODUCTION 
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0^ i Charged polymers have drawn a considerable amount of attention in the past years both theoretically and experi- 
mentally JPH| • This is due, on one hand, to their wide range of industrial applications in processes involving charged 
colloids mM and on the other hand, to their resemblance to water soluble bio-polymers such as proteins and DNA. A 
Q distinction is made between polyelectrolytes, with all charges having the same sign, and polyampholytes with positively 
as well as negatively charged monomers |6|-|lc||. For the former, the electrostatic interactions are repulsive and long 
ranged. For the latter, the repulsion of like-charges competes with the attraction of opposite charges, resulting in a 
q ' complex behavior which depends on the net charge of the chain. 

In spite of extensive investigations of polyelectrolytes including the pioneering works of Katchalsky et al. and 
Manning |fT2| , polyelectrolytes are much less understood than neutral polymers. For example, there exists a debate 
regarding the persistence length and the chain conformations for single chains in dilute solutions p^-jle| . This is due 
to the delicate interplay between the chain connectivity and the long range nature of electrostatic interactions. 
Semi-dilute polyelectrolyte solutions where the chains interact with one another have been studied as well 
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6|,fl 7 19 2fJ]. In the so-called "blob" picture, scaling laws are derived by singling out the most dominant inter- 



action at different length scales. Another technique is the Random Phase Approximation (RPA) [Ell E6| used for the 



calculation of the structure factor S(q) as can be measured by scattering experiments 27|26 

In addition to bulk properties, special attention was directed to adsorption experiments of polyelectrolytes onto a 
single charged surface p9|-|32|] , and between two charged surfac es [^3|-^6| . On the theoretical side, discrete models 
have been employed for which the chains are placed on a lattice p7pH0|7This approach has been used by Bohmer et 
al. [JO] to calculate force curves between two charged surfaces containing a polyelectrolyte solution. 

Another approach is a continuum one [^2|-^6| where the charge densities, monomer densities and electric field are 
treated as continuous functions of the local position. Varoqui et al. Jl3|,fl4j investigated polyelectrolyte adsorption 
onto one surface, while Podgornik §47[ calculated inter-surface forces. In recent works |4q] non-linear excluded volume 
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interactions have been included and scaling laws characterizing the adsorption of polyelectrolytes are proposed (within 
mean- field approximation). 

In the present work the continuum approach is extended to random (heterogeneous) polyelectrolytes. We study 
several statistical charge distributions (i.e. annealed, quenched and permuted) corresponding to different physical 
situations. Polyampholytes with positive and negative charges are studied as well. One of the goals of the present 
work is to take explicitly into account several characteristics of polyelectrolytes such as the connectivity of the 
polymer chains, the non-electrostatic monomer-monomer interactions and the Coulomb interactions between charged 
monomers, counter-ions and co-ions. 

The paper is organized as follows: in the next section we present a general formalism for calculating the free energy 
of randomly charged chains. This formalism is applied in Sec. Ill to derive mean-field equations for the various charge 
distributions (including polyampholytes). The reader who is not interested in the technical details can skip these two 
sections and go directly to Sec. IV where the mean-field results are summarized and the various charge distributions 
are compared. The random phase approximation (RPA) is used in Sec. V to calculate the structure factor S(q) of 
polyelectrolytes and polyampholytes in good and bad solvents. 



II. GENERAL FORMALISM 



Let us consider a semi-dilute solution of polyelectrolytes in a good solvent in presence of salt (electrolyte). The 
system is schematically drawn on Fig. 1. In the model, based on the Edwards' formalism [Q, the microscopic degrees 
of freedom are the monomer positions {r;(s)}, where s € [0,iV] is a continuous index along a chain of N monomers 
and I = 1, M is the label of the M chains in solution. The positions of the small co-ions (counter-ions) are denoted 
by R+ (Rj) where i — 1, N + and j = 1, N~ are, respectively, the indices of the monovalent positive (+e) and 
negative (— e) ions (see Fig. 1). For simplicity only a symmetric 1:1 electrolyte is considered. 

The partition function Z is then expressed as a path integral over all possible configurations: 



Z = J Vnis) VR+ VR- exp jTda rf(s)j 

x exp f-^/? J dr dr' p c (r)v c (r - r')p c (r') -^v J dr pi 



(1) 



where J T>g denotes the functional integral over the function g(r). The first term in the exponent is the Wiener 
measure representing the connectivity of the polymer random walk ]49|], where r;(s) is the derivative of r;(s) with 
respect to the monomer index s and the Kuhn length a is the effective monomer length. In writing this term, we have 
assumed that the charged chains can be modeled as flexible chains. This assumption can be more easily justified for 
semi-dilute solutions of weakly charged chains. 

The second term is the electrostatic interaction term where /3 = 1/ksT is the inverse thermal energy, v c (r) — l/er| r| 
is the Coulomb interaction, e is the dielectric constant of the solution and /5 c (r) is the local charge density operator 
including all charges in the solution (charged monomers, co-ions and counter-ions): 

M N N+ N~ 

p c (r) = J2 dsq l (s)6(r-r l (s)) + J2eS(r-Rt)-J2 e6 ( r - R j) ( 2 ) 

1=1 •'° z=l 3=1 

where qi (s) is the random variable denoting the charge carried by the s monomer along the I chain. In the next section 
we will consider several charge distributions for qi(s). For example, in the smeared model, qi(s) = pe is a constant 
independent of the position s along the chain and the chain index I. 

The last term in eq. [j] is the excluded volume repulsion between monomers, where v ~ a 3 is the excluded volume 
parameter and 

M <>N 

ft»W=E/ d S <5(r-rKs)) (3) 

is the local monomer concentration operator. 

It is possible to integrate out the microscopic degrees of freedom in the partition function eq. [I] by introducing two 
pairs of collective coordinates: (i) the local monomer concentration p m (r) and its conjugate field ^5 m (r); and, (ii) the 
local charge density p c ( r ) and its conjugate field y c (r). This is done using the following identity 
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1 = J Vp m 5(p m - p m ) = J Vp m Vtp m exp ^iv J dr p m v? m -w^^J ds <p m [*l(s)]\ 



(4) 



for /f m (r) and an analogous identity for p c ( r )- We will see below that i<p c ( r ) is in fact the electrostatic potential. The 
functional integrals over p m (r) and /? c (r) are Gaussian integrals and are evaluated exactly, leading to 

V^ m T>ip c exp (- J dr j/?^|V<^ c | 2 + ^v&Vj C+ C- C P (5) 

where (± = (J dr e Tl/3ev °( r )) A ' ± and £ p are, respectively, the partition functions of the (small) co-ions, counter-ions 
and polymer in the presence of the external fields i<p c ( r ) and i</3 m (r). In the thermodynamic limit, where and the 
volume V become large while the bulk concentrations = N^/V remain fixed, £± become (up to a normalization 
constant): 



(6) 



1 + 1 Jdr (e* i0ev °W - l) " _ » exp^c± ^dr {e^^'M - l}^ 
The partition function of the polymer chains £ p in the presence of the two external fields is: 

F rf (s) + ivip m [n{s)] + i/3qi{s)ip c [ri(s)]^j (7) 

Note that the calculation of ( p depends on the specific charge distribution for qi(s). 



C P = Jvr l{s) exp^g^d, {A r j 



III. MEAN-FIELD EQUATIONS 

We apply the general formalism introduced above to different monomer charge distributions: smeared, annealed, 
permuted, quenched and polyampholytes. These charge distributions are applicable to different experimental systems. 

A. "Smeared" Polyelectrolytes 

The simplest model of a charge distribution is the so-called smeared polyelectrolyte. For a polyelectrolyte with a 
fraction p of its monomers being charged, this model assumes that each monomer carries a uniform fractional charge 
pe, where e is the electron charge. Namely, qi(s) — pe for any monomer s on any chain Without loss of generality 
we shall assume that the polymer is positively charged. 

It is possible here to use a well-known analogy from quantum mechanics to calculate the path integral of eq. M. 
The partition function is analogous to a Feynman integral of the Hamiltonian TL — — + iv(p m + i[3pe(p c with 
imaginary time t — > is. Thus, each eigenstate <f> v contributes a term of the form exp(—MNE v ) where MN is the 
total number of monomers in the solution. In the thermodynamic limit N — * oo, the ground state dominates over all 
other eigenstates |5C]], and Cp reduces to: 

C P « e~ MNE ° = exp (- J dr |^|V0| 2 + l ^ m (r)0 2 (r)+z/?pe^ c (r)0 2 (r)- Alp 2 (r)|j (8) 

where Eq is the ground-state energy, <j)(r) is the renormalized ground-state eigenfunction and /i p is a Lagrange 
multiplier added in order to ensure the normalization of the wave- function 2 = y J dr </> 2 (r), <^> 2 being the bulk 
monomer concentration. The polymer analog of the wave-function <f>(v) is usually referred to as the polymer order 
parameter, and the local monomer concentration per unit volume can be shown to be p m (r) = 2 (r) 

The field tp m can now be integrated out leaving a functional integral only over ip = itp c , and the partition function 
reduces to: 

Z = JviJj exp(-/3F s ) (9) 

where 
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0F S = J dr |-/?^-|VV| 2 + ^c±(l-c 



=F/3ei/>(r) 



+ y|V^| 2 + ^0 4 (r) - Mp 2 (r) + /3pe</>(r)</> 2 (r) j (10) 

Note that i/>(r) is identified as the electrostatic potential. 

Within a mean-field approximation, the functional integral is dominated by the saddle point given by the condition 
5F/5ip — 0. This results in a Poisson-Boltzmann (PB) like equation for the electrostatic potential, which includes all 
charge sources 

VV(r) = — c b sinh(/3e^) - — (W> 2 - p&eP**) (11) 
e e 

The right hand side represents the local charge density. The first term is the symmetric 1:1 electrolyte contribution, 
the second term is associated with the (positive) charges on the polymer chains, while the last term comes form the 
counter- ions which dissociated from the chains. Note that charge neutrality implies — * Cb and — > Cb + P</> 2 - 

Since <^(r) is the ground state eigenfunction, it satisfies the variational equation 8F/5<p — 0, yielding a self-consistent 
field (SCF) equation for the polymer order parameter: 

^V 2 </>(r) = «(0 3 - 4>U) + Ppe^ (12) 
6 

This is an equation for the density of polymer chains in an external electric potential tp and with excluded volume 
interactions characterized by the parameter v. In the above equation n p has been substituted in order to have the 
correct bulk limit: ip — > and <fi 2 — > (j)^. 



B. Annealed Polyelectrolytes 

The derivation presented above for the smeared case can be extended to more realistic annealed charge distributions 
where charges can dynamically associate and dissociate from the chains. The annealed model describes an experimental 
system where the monomers have weak acidic (or basic) groups. The pH of the solution controls the degree of 
association/dissociation of ions on the polymer chain. 

Assuming no charge correlations along the chain, the monomer charge distribution f[qi(s)] is defined as the nominal 
probability of the monomer s along the chain I to have a charge qi(s). 

f[m(s)]=J2Pj s (<n( s )- z 3 e ) ( 13 ) 

3 

where the randomly charged chain is described by a set of valencies {zj} (zj = 0, 1, 2, . . .) with normalized probabilities 
{pj} such that ^2jPj = 1. Here, we concentrate on a simple example of a polyelectrolyte for which each monomer 
can be either positively charged (z = 1 and q = +e) with probability p or neutral with probability 1 — p. The charge 
distribution for each monomer is 

/[©(*)] =pS(qi(s)-e) + (l-p)6( qi (s)) (14) 

It is important to note that for annealed polyelectrolytes, the partition function has to be averaged with respect to the 
monomer charges since they are in thermal equilibrium with the reservoir. As a result the electrostatic contribution 
of one monomer in ^ p (eq. |?]) becomes: 

^ e -*/3'3i(s) l Pc[r;(s)]^ _ j _ p _|_ pg-*/3eip c [rj(a)] 

where the average is taken over the charge distribution, (0) p = J 0(q)f(q) dq, and the annealed free energy becomes 
/3F a = J dr j-^IV^ + ^c^l-eWW) 

+ y|V0| 2 + ^ 4 (r)- Mp ^(r)-^(r)log(l-p + pe-^«)j (16) 
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The variation of eq. with respect to ^(r), leads to a PB like equation 



where 



VV(r) = ^ Cb sinh(/3e^) - ^ (p a tf - p<&^) 



(17) 



(18) 



can be interpreted as the annealed charge probability and depends on the local electrostatic potential. Similarly, the 
variation of eq. |l6| with respect to <p(r) gives a modified SCF equation: 



(r) = v( 



>l0g(l -J3 a +PaC 



(19) 



Equations [I?] and [l^ are the annealed equations for the charge distribution eq. |l4], similar to eqs. 11 and 12 for the 
smeared case. 



C. "Permuted" Polyelectrolytes 



Another variant of the annealed case is the permuted model for which a fixed number of charges pN are free to 
move along each chain without dissociating from it. Thus, the total charge on the chain remains constant. This is 
introduced into the annealed model by adding a constraint in eq. [?] in order to keep the total charge on each of the 
chain fixed: 



N \ , f N 

ds qi(s) — Npe = / dui exp iui / ds [qi(s) — pe] 



(20) 



where ui is a Lagrange multiplier related to the fixed charge constraint of the chain I, Assuming ground state 
dominance, the polymer partition function becomes: 



du exp 



dr 
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|V(/>| 2 + ivtp m (r)(j> 2 (r) + ipeu4> 2 (r) 



Mp^(r) - 2 (r) log (l + e-^M+- 



(21) 



where the index I is dropped from the functional integral over ui since the constraint is satisfied separately on each 
chain. Note that the fraction p of charged monomers is introduced through the constraint (eq. |20| ). In order to carry 
out the integration over u, we use the identity (similar to eq. 0) 



1 = J T>g(r) %(r) - u) = J Vg(r) Vh{v) exp (i J dr h(r)g(r) - i J dr 



h(r)u 



(22) 



A sequence of saddle point approximations for g(r), h(r) and u leads to the following mean-field free energy for the 
permuted case: 



0F p = pF s + J dr {p\ogp + (l-p)lo g (l-p)}(£ 2 (r) 



(23) 



This correction term represents the translational entropy of the charges on the chain p3]. However, the last term 
being quadratic in (j) and independent of ip, it only shifts the chemical potential /i p without affecting the differential 
equations [ll] and [l^. 
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D. Quenched Polyelectrolytes 



In the quenched model the charge distribution is frozen. Experimentally, this corresponds to heterogeneous copoly- 
mers with a random sequence of charged and neutral monomers. The specific sequence of each copolymer is determined 
during the polymerization stage and represents one possible realization of the random distribution. Various physical 
quantities are then calculated by averaging over this random distribution. Instead of averaging the partition function, 
one should average the free energy itself over the random charge distribution. A standard method in quenched systems 
is the replica trick [ p2| based on the following identity: 

(logZ) = lim (24) 

where (. . .) indicates an average over the quenched disorder and the n-th power of the partition function introduces n 
replicas of the system all having the same monomer charge distribution. Generalizing the above approach we obtain: 

(Z n ) = Iv^Vtf expl-V /dr {[3^-\X7tf\ 2 + V c±(l - e^°) + l v (<) 2 1 ] x (Cn) (25) 



V^Vy* exp(-^ J dr |^|V^| 2 + ^c±(l 
\ a=l y ± 



where 



',<,.) = J Vv?{s) expf-^^jf ds [^(rf(s)) 2 W<v;;,[ r ;'(,s)] 

/ M r N \ 

xexpf^y ds log^l-p+pe-^Ea^l'rWlj J (26) 

and the superscript a = l,...,n is the label of the replica a. A mean-field estimate of can be obtained by 
using again the quantum-mechanical analogy. An additional complication is that the effective Hamiltonian here is a 
many-body one as it couples different replicas: 

n ( 2 ~i 

E \-\ V « + w < - log (l - P + pc-^« (27) 



Hi. 

Assuming no replica symmetry breaking, we use a Hartree approximation to express the ground state of the many-body 
eigenfunction <f) n ({r a }) = []™ =1 </>(r Q ), in terms of single-body eigenfunctions 4>(r a ): 

C„ « e- ME » = exp f-nM ^ dr |y |V^| 2 + iv^ m 4> 2 - M P <^}) 

x exp y dr Q f[ <j> 2 {v a ) log(l - p + pe" i/3e £a ^(O^ (28) 

The path integrals over <Pm( r ) can t> e evaluated exactly, whereas the path integrals over cp" (r) can be approximated 
by their saddle point values. Since we assumed no replica symmetry breaking, all saddle point functions are identical: 
icp"(r) = ip(r). When we take into account the coupling with the polyelectrolyte reservoir, this equation reduces to 
the annealed PB equation (eq. [I?]). Similarly, the SCF equation reduces to eq. [H]. 



E. Annealed Polyampholytes 

Finally, the mean-field formalism can be generalized in a straightforward way to treat polyampholytes, consisting 
of negatively and positively charged monomers. A general polyampholyte is described by a set of valencies {zj} 
(zj = 0, ±1, ±2, . . .) with probabilities {pj}, where J2jPj = 1- The statistical distribution of the charges can be 
either annealed, quenched or permuted. For simplicity, we consider only the annealed case where the contribution of 
the charged monomers to the free energy becomes 

= 



J dr 4> 2 {v) log ^p e- pz ^ (r) (29) 



G 



and should be compared to the last term in eq. [T(]. The modified PB equation (eq. |T^) is now: 



V 2 ^(r) = ^c b sinh(/3e^) - ^ {z a <j? - zcge^) (30) 
where z = ^2jPjZj is the average monomer charge (in units of e) and 

is the annealed (weighted) monomer valence, which is nothing but the Boltzmann average of the monomer charge 
distribution. Finally, the SCF equation (eq. |l^) is: 

yV 2 ^(r) = - fy) - <j>\og [5>e-^ J (32) 



IV. DISCUSSION OF THE CHARGE MODELS 



In the previous section we derived the free energies and mean-field equations for different charge distributions, by 
path integral methods. 

The smeared case: The smeared free energy can be separated into uncharged polymer and Coulombic contributions 
g|, F s = F pol + F ol . The polymer part is: 



F pol = k B T J dr | ^|V0| 2 + \v^{v) - Mp 2 (r) j (33) 

This is the Edwards free energy for polymer solutions expressed in terms of the polymer order parameter </>(r), which 
is the square root of the monomer concentration, p m = cj) 2 . The first term represents the stiffness of the polymer 
chains, the second takes into account the excluded volume in good solvent conditions, while the last term represents 
the coupling to a reservoir of monomers with chemical potential /i p . 
The Coulombic free energy reads: 



Fei = J dr l-l-\Vi>\ 2 + k B Tj2ct(l- 



e TWW) + pe^( r )^(r) (34) 



This free energy contains the electrostatic interactions between all charges (small ions and charged monomers), as 
well as the translational entropy of the small ions in solution. In the above expression the independent fields are the 
electric potential ip and polymer order parameter <j>. The small ion concentrations are uniquely determined by the 
electric potential through the Boltzmann weight: c^iv) = exp(=F/3e-0(r)). 

The same free energy F e \ can be obtained from a more direct approach as will be shown for the smeared case. It is 
convenient to express the free energy in terms of the total charge density: p c = ec + — ec~ + peip 2 

# cl = Ifdrdr' 



2 1 elr- 



+ k B T J dr {c+(logc+ - l) + c - (logc- -1)}- J dr { A i+c+ + ft_c - } (35) 

The free energy _F c i of eq. |34| should be distinguished from F c i. The former depends on ip and <fi while the latter 
depends on and <p. For the latter, the electric potential can be defined as 



and it satisfies the Poisson equation 
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(37) 



Minimizing F e \ with respect to c we obtain the equilibrium charge distribution of the small ions: 

c ± (r) = c±exp(T/3e^(r) (38) 

where fi± = fc^riogc^. 

Substituting the above equilibrium condition back into eq. ^ we obtain 

F ei = ^ J dr pecp 2 ip - i / dr (ec+ - ec~)V> 

-k B T J dr {c+(e-^ - 1) + c^e^ - 1)} (39) 

From the Poisson equation ^ we can express c + — c~ in terms of V 2, and </>. Integration by parts of the term ipX7 2 il> 
yields exactly the first term of eq. |34| with the correct negative sign |^5|,^l| . 

The annealed case: The second type of charge distribution is the annealed one where each monomer can be either 
charged or neutral with bare probabilities p and 1 — p, respectively. Its free energy eq. |l6| is similar to the smeared one 
except for the coupling term between the charged monomers and the local electric potential. The difference can be 
understood in the following way: in the annealed case different charge configurations will contribute to the free energy 
and one needs to trace over those configurations in the partition function before the free energy is calculated. In the 
smeared case there is only one charge configuration where every monomer is assigned a fractional charge pe, whereas 
in the annealed case p represents the bare probability of dissociation and p a — pexp(— f3eip)/(l — p + pexp(— (3eip)) is 



the effective probability as can be seen in eqs. 17 and h 



In experiments the effective probability p a is related to the pH of the solution via 

P H = P K + log 10 (40) 

^ Pa 

where pK = — log 10 K and K is the dissociation constant. For example, in the case of weak alkaline monomers 

AOH^A+ + OH" (41) 

Ko is given by 

[A+HOH-] 

Ko - [AOH] (42) 

From eq. ^ one can easily obtain p a as a function of the pH. 

At low electrostatic potentials \{3eip\ <C 1, the annealed free energy can be expanded in powers of ip. The first term 
is equal to the smeared free energy F s , while the next term is always negative 

/?F a ~/?F s -^(l-#V|dr^WW < PF S (43) 

The fact that F a < F s is related to the convexity of the free energy. Indeed, the annealed case has more degrees of 
freedom and allows a better minimization. 

The permuted case: The permuted model is a variant of the annealed case. It models either mobile charges which 
can hop along the chain or charges which associate and dissociate while keeping the total amount of charge fixed on 
each chain. 

The free energy for the permuted case can be written as the smeared free energy eq. [l^ plus an additional term of 
entropic origin: 

Fp = F s + ksT J At {p\ogp + (1 - p) log(l - p)}4> 2 {r) (44) 

The correction is due to the translational entropy of the charges along the chains. It amounts to a shift in the polymer 
chemical potential and thus does not affect the mean field equation. 
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The correspondence between the permuted and the smeared models was not emphasized in previous works. It 
can be interpreted as a tendency of the charges in the permuted model to spread uniformly along the linear chain. 
However, as was discussed earlier [|3|, changing the pH of the solution (e.g., by titration) can lead to non trivial 
dependence of \x p on the physical parameters, since the reservoir concentration changes in a titration process and will 
affect /i p . 

The quenched case: In order to obtain the equilibrium state of chains with frozen (quenched) charge distributions, 
the free energy has to be averaged over all possible charge configurations. As was shown in the previous section, 
annealed and quenched polyelectrolytes in contact with an infinite reservoir of chains have the same mean field free 
energy. 

The physical meaning of this result can be explained in the following way (53): when quenched polymers are 
allowed to exchange with a bulk reservoir, containing all possible configurations, the system picks up the optimal 
configurations from the bulk. When the polymers are not coupled to an infinite reservoir the annealed and quenched 
cases arc different. Note that the dynamics of annealed and quenched polyelectrolytes can differ considerably, but 
this is beyond the scope of the current work. 

The annealed polyampholyte case: This situation corresponds to monomers which can carry a positive or 
negative charges of valency Zj = 0, ±1,±2, ... with probability pj. As was derived in the previous section the 
electrostatic part of the free energy is given by: 



pel 

pa 



= -k B T J dr </> 2 (r)log ( ^pp-?*!^ J (45) 



For low electrostatic potentials, an expansion of the above polyampholyte free energy yields: 

1 
2 



F£ ~ ze j dr ip(r)(j) 2 (r) - -pa 2 e 2 / dr ip 2 {r)<f> 2 (r) (46) 



The first term is linear in ip and reduces to the smeared contribution where ze replaces pe as the average monomer 
charge. The second term is a negative correction (as expected) which depends on the statistical variance of the charge 
distribution 



J2 _ „ ~2 



X>*j) (47) 



/ J 1 

3 3 

This term affects the monomer- monomer correlations and S(q) as discussed in the following section. 
The general distribution treated here has a few simple and useful limits: 

1. The smeared model (Sec. III. A) is exactly recovered for px = 1 and z\ = p yielding z = p and a 2 = 0. 

2. The annealed model (Sec. III.B) is obtained for p% = p, z\ = 1 and pi = 1 — p, z-i = 0; z being simply p, z a 
being p a , and the variance a 2 = p(l — p). 

3. A trimodal charge distribution can be a good representation for some polyampholyte systems. 

f(q) = p+6(q - e) + p S(q) + pS(q + e) (48) 

where po = 1 — p+ — p—. Each monomer can be either positively charged, negatively charged or neutral 
with probabilities p + , p_ and po, respectively. The trimodal distribution is characterized by two independent 
parameters. These can be the first two moments: the average z = p + — P- and the variance ("width") 
a 2 = p_|_(l— p + )+p_(l— p_) + 2p + p_. Increasing the weight of po means that the polyampholyte becomes weakly 
charged, while increasing the weight of p+ (or p_) means that the polyampholyte becomes more asymmetric 
and resembles more a true polyelectrolyte. 

4. A bimodal distribution of polyampholytcs where each monomer carries either a +e charge or — e one, with 
probabilities p and 1 — p, respectively. 

f(q)=pS(q-e) + (l-p)S(q + e) (49) 

with z — 2p — 1 and a 2 = 4p(l — p). Note that for this fully charged polyampholyte the same parameter p 
characterizes both the average z and the variance a 2 , so that a 2 = 1 — z 2 . The distribution is symmetric around 
p = 1/2. As |p — 1/2| increases, \z\ increases, a 2 decreases, and the polyampholyte resembles more and more a 
polyelectrolyte with a net charge. 
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It is of interest to consider explicitly the symmetric bimodal case mentioned above having no net charge p = 1/2, 
z — and a 2 — 1 is maximal. The PB equation (eq. ^0|) can be written as 

o7TP 47TP 

V V(r) = c b sinh(j3eip) + <j> 2 tanh(^eV>) (50) 

e e 

and the SCF equation (eq. ||^) is: 

^V 2 ^(r) = «(0 3 - 4>&) - ^log [cosh(/3e^)] (51) 

The second term in eq. |5^ represents the contribution of the charged monomers to the local charge density. At low 
potentials \(3eip\ *C 1 the polymer charge density is e<fi 2 tanh(/3e - 0) ~ e<f> 2 smh((3eip). Comparing, in this limit, the two 
terms on the right hand side of eq. |5^, the polyampholyte can be viewed as a symmetric electrolyte || whose bulk 
concentration is not a constant but determined by the local monomer concentration. For a monovalent electrolyte, 
the local concentrations of positive and negative ions obey a Boltzmann distribution c (r) = Cb exp(=p/3ei/'), where 
the equilibrium distribution is achieved by exchange of ions with the reservoir. For the annealed polyampholytes, 
the concentrations of positive and negative monomers j0j^(r) behave in a similar way (for weak potentials), Pm( r ) — 
^(j) 2 exp(=p/3e , (/'), but the mechanism is different. For the latter, the equilibrium distribution is achieved by ionizing 
the monomers, and this process is limited by the monomer concentration, whereas for the former case, the reservoir 
contains an infinite amount of ions. This difference becomes evident at high potentials where tanh(/3e^>) saturates to 
+1 or —1, depending on the sign of the potential. Under these extreme conditions the polyampholytes are no longer 
neutral. Instead, they are fully ionized. 



V. STRUCTURE FACTOR WITHIN RPA 

Density-density correlations are measured in scattering experiments ^?],^8| and can be calculated using the random 
phase approximation (RPA) [pl|— p6| . This is done by considering small fluctuations of the homogeneous bulk state. 
Since we are interested in the monomer- monomer density correlations, we do not perform the integration over p m (v) 

in eq. |[ and express the partition function in terms of the three fields p m (r) — <fi 2 + Sp m (r), (p m (r) = tp^ + 5<p m (r) 

and </? c (r) = <pc + S(p c (r), As introduced in Sec. II, ip m is the conjugate field of the local monomer concentration p m 

and ip c is the electric potential conjugate to the charge density p c . Note that ipi°^ — in the bulk. 

The free energy is then expanded to second order in these fluctuations. Since the linear terms in Sp m , <5iy9 m and 5ip c 
cancel out, the first non-zero corrections are of second order. 

It is more convenient to write this expansion in Fourier space: 

^ m (r)=y^3 5p m (q)e iq ' r (52) 

and similarly for 5ip m and 5ip c . 

A. Smeared 5(g) 
Using the smeared free energy, the expansion yields 

(3F s {8p m ,5ip m ,5(p c ) ~ /t^t { vSp m (-q)Sp m (q) + v 2 (/)lSo(q)6tp m (-q)Sip nl (q) 




+P^-[q 2 + k 2 s + pK 2 p S (q)\ 5(p c (-q)5(p c (q) 
+Ppev4>lS Q {q)6(p m (-q)Sip c (q) - ivSp m (-q)S(p m (q)\ (53) 

where the integral is over the wavevector q, and So(q) is the structure factor of Gaussian chains. For chains of length 
N, it is equal to S (q) = ND(\a 2 q 2 N) where D(Q = (2/( 2 )(e" c + £ - 1) is the Debye function g. For infinitely 
long chains N — > oo, the structure factor is independent of the chain length, 5 , ^~ 1 (g) = j^a 2 q 2 . 
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The Debye-Hiickel screening length has two contributions: 

K tot = K l+ P K l = 87r/ B c b + Airl B p4>b ( 54 ) 

the first contribution arises from the 1:1 symmetric electrolyte while the second one comes from the polymer counte- 
rions, Ib = e 2 /efc^T being the Bjerrum length, equal to about 7A at room temperature. 

The Fourier transform of the monomer-monomer correlations (Sp m (r)6p m (0)) is proportional to the experimentally 
measured structure factor S(q). It can now be calculated as a Gaussian integral giving: 

S s\q) = 77 7^77 TT7 = S 0~\<l) + v $o + 9 ~ 9 P 7 ( 55 ) 

(<5p m (-q)<5p m (q)) Yh q 2 + K 2 s + P Kl 

Similar expressions were obtained in previous works [pl|-p5[, using somewhat different derivations. 



B. Annealed S(q) for Polyelectrolytes and Polyampholytes 

Repeating the above calculation for the annealed case leads to the same form for S(q) where n 2 ot is replaced by 
K tot +P(1 - P) K l in eq. H- 

S,\a) = S;%) + v<& + q , + /lf {2 _ pWp (56) 

The enhanced screening here is due to the additional annealed degrees of freedom of the charges on the polymer 
chains. Local fluctuations of the monomer charge density effectively increase the local ion concentration leading to 
stronger screening. 

The structure factor is also calculated for annealed polyampholytes. It amounts to replacing n 2 at byK 2 + (|z|+cr 2 )K 2 
and p 2 n 2 by z 2 k 2 in eq. |55| leading to: 

= Vfe) + ^ + (f + f ,^ l + a . H (57) 

For neutral (symmetric) polyampholytes, z = 0, expression^ is the same as S~ 1 (q) of neutral polymers. We note 
that this is an outcome of the RPA which neglects higher order charge correlations. However, the mean field equations 
themselves as well as the free energy ^5] depend on the charge distribution and, in particular, on the variance a 2 . 

At high charge fraction \z\ and low salt concentration, the structure factor S(q) exhibits a peak at a finite wavenum- 
ber qo > satisfying 

[ql + k 2 s + (\z\ + a 2 )K 2 ] 2 = 12z 2 K 2 p /a 2 (58) 

This peak is characteristic of polyelectrolyte solutions at low salt concentration. Since S~ 1 (q) is the energy of the 
q mode density fluctuations, a maximum in S( q) corresponds to the lowest energy fluctuation. Here qo > results 
from the competition between the first term in |37|, originating from the polymer elasticity, and the last term due to 
the electrostatic interactions screened by the small ions. For q values K < q < a -1 , the peak qo can be observed. For 
example, in the absence of salt, the condition for having a peak at qo > is 

<f!<^-l (59) 
\z\ K p a 

The upper bound is bigger than zero as long as k" 1 > a recalling that k~ 1 = 1/ y/^rds^ and a is the monomer 
size. The above inequality is satisfied for highly asymmetric polyampholytes, resembling polyelectrolytes (see case 3 
of the preceding section). In the opposite limit, low \z\ and/or high a, the random polyampholyte behaves essentially 
as a neutral polymer (go — * 0). For fully charged polyampholytes (case 4 of the preceding section), a 2 /\z\ = (1 — 
z 2 )/\z\ yielding the same conclusion. As long as the net charge \z\ is large enough, the polyampholytes resemble 
polyelectrolytes and S{q) will exhibit a peak at finite qo- 

In Fig. 2 the structure factor S(q) is plotted as function of the wavenumber q for smeared polyelectrolytes (a 2 = 0) 
for different salt concentrations and charge fractions z. The structure factor at wavenumber q = is proportional to 
the osmotic compressibility. As depicted in Fig. 2, S(q = 0) increases upon addition of salt, while the peak position at 
qo shifts to smaller wavenumbers until the peak disappears and S(q) becomes a monotonous decreasing function. The 
inset of Fig. 2 shows the effect of decreasing the average monomer charge \z\. The peak increases and shifts towards 
smaller values of q. 
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C. Mesophases in Bad Solvent 



We end this section by examining polyelectrolytes and polyampholytcs in bad solvent conditions. The excluded 
volume parameter v is negative leading to collapse (and segregation) of the chains, and higher order virial terms have 
to be included in the free energy. Assuming that the third order virial coefficient w is positive, with a contribution of 
\wp^ to the free energy, we obtain 

= S^ iq) - Wl + W <f>i + — (60) 

An analogous expression was obtained for smeared polyelectrolytes j2l] , p2|j25f and is generalized here to annealed 
polyampholytes. 

An instability of the homogeneous phase is determined by 5 ,_1 (qo) = 0, indicating a mesophase separation, where the 
size of the micro-domains is Ao = 27r / ' q^. In Fig. 3, the line marking the instability of the homogeneous (disordered) 
phase is plotted for different salt concentrations. At high ionic strength dilute polyelectrolyte solutions becomes 
unstable. This reflects the macroscopic phase separation of neutral (screened) polymers in a bad solvent (qo — > as 
is seen in Fig. 2). 

In Fig. 4 we show the dependence of the micro-domain region on the variance of the charge distribution, a 2 . 
One important consequence of eq. [6^ is the equivalence of electrostatic screening induced by fluctuations in the 
polyampholyte charges and screening by added salt with effective bulk concentration c£ ff = \d 2 ^\- For example, for 
4>l = 1(T 6 A- 3 , cf/a 2 = 0.8mM. This effect is stronger in polyampholytes than in annealed polyelectrolytes, since 
in almost symmetric polyampholytes the excess charge is small, \z\ <C 1, while the variance of the charges a 2 can be 
close to one. 

It is instructive to look at the instability of the homogeneous system towards macro-phase separation as a function 
of the polyampholyte average net charge. We emphasize that in some cases this macro-phase separation occurring at 
q = is preempted by a mesophase at qo > 0. Nevertheless, let us consider the change in the second virial coefficient 
v ~ T — Tg due to the electrostatic interactions where Tg is the ^-temperature in the absence of electrostatic 
interactions. We note that for z = there is no electrostatic contribution to S(q) in our RPA calculation. For 0, 
the q — instability will occur for 

M« " «f + (wt°K (61) 

For fully charged polyampholytes (taking the bimodal distribution: Zj — ±1 and \z\ < 1), the 0-temperature in the 
presence of charged monomers is 

2cb/0b + \ z \ + 1 - z 

This behavior is similar to the one found for single chains . In Fig. 5 we plot the dependence of the 0-temperature 
on the net charge z for different salt concentrations. As expected, the (9-temperature is a decreasing function of the 
net charge \z\. As can be seen from Fig. 5, addition of salt extends the bad solvent regime to higher values of \z\. At 
high salt concentrations, the electrostatic interactions are screened and Tg(z) is very close to its pure value Tg. 



VI. CONCLUSIONS 



In this work we have studied bulk properties of charged polymers in aqueous solutions in the presence of added 
salt. Starting from a path integral formalism which takes into account the chain connectivity, short range and electro- 
static interactions we derived mean-field equations describing the behavior of polyelectrolytes and polyampholytes in 
solution. We compared several models for the statistical charge distribution corresponding to different experimental 
realizations. The simplest and most frequently used model is the smeared one where charges are uniformly distributed. 
The permuted model, where the charges are mobile along the chain, is found to be equivalent to the smeared model 
except for a constant shift in the monomer chemical potential. This shift has to be taken into account in titration 
experiments. 

The annealed model was found to have a lower free energy than the smeared one. This is related to the additional 
degrees of freedom. The effective dissociation depends on the local electric potential. At thermodynamic equilibrium, 
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the quenched case is found to be equivalent to the annealed one, as long as the system is in contact with an infinite 
reservoir (bulk) of polyelectrolyte chains. 

Annealed polyampholytes are characterized by their net charge and variance. We find different behavior for symmet- 
ric polyampholytcs (no net charge) as compared to asymmetric ones (closer to polyelectrolytes). At low electrostatic 
potentials, all the above mentioned models have the same limiting behavior. 

The monomer-monomer structure factor S(q) is calculated within the Random Phase Approximation for annealed 
polyelectrolytes and polyampholytcs. The electrostatic screening depends not only on the salt and counterions but 
also on the variance of the annealed charge distribution. 

A peak in S(q) for polyampholytes at finite wavenumber is shown to appear at a high net charge and/or low variance 
indicating polyelectrolyte-like behavior. For bad solvent conditions, the variance enhances the tendency of the system 
to undergo a mesophase separation. 

Finally, it is worth mentioning that the present study can be further extended to treat polyelectrolytes and polyam- 
pholytes in restricted geometries and close to charged surfaces |u],[l3 45- 4^]. In particular, it will be interesting to 



address the question of how the different charge distribution is coupled with the polymer adsorption onto a single 
surface and the forces exerted by polyelectrolytes and polyampholytes between planar, cylinders or spherical surfaces. 
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FIGURE CAPTIONS 



Fig. 1: Schematic view of a polyelectrolyte solution. The monomer coordinates are r;(s) where I = 1, ...,M labels 
the polymer chain and s £ [0, N] is a continuous index along the chain. The small ion coordinates are and 
R~ where i = 1, ...,N + and j = 1, ...,N~. 

Fig. 2: The effect of salt concentration on the structure factor S(q) of polyelectrolytes, eq. |57|. The parameters used 
are: 4> 2 = 10 _6 A -3 , a = 5A, v = 0.1a 3 , z = 0.1 and a 2 — 0. The salt concentration is: Cb = ImM (solid curve); 
Cb = 2mM (dots); Cb = 5mM (short dashes) and Cb = lOmM (long dashes). The inset shows the effect of the 
average monomer charge at small salt concentration Cb = O.lmM. The different curves correspond to: z = 0.10 
(solid curve); z — 0.03 (dots) and z = 0.01 (dashes). 

Fig. 3: The effect of salt concentration on the stability line of the homogeneous phase. The parameters used are: 
a = 5A, w = 5 x 10 5 A 6 , z — 0.01 and a 2 = 0. The salt concentration is: Cb = (solid curve); Cb = ImM (short 
dashes); Cb = 1.5mM (dots and short dashes); Cb = 2.5mM (long dashes) and Cb = 5mM (dots and long dashes). 

Fig. 4: The effect of the variance a 2 on the stability line of the homogeneous phase. The parameters used are: 
a = 5A, w — 5 x 10 5 A 6 , z — 0.01 and Cb = O.OlmM. The different curves correspond to a 2 = (solid curve); 
a 2 = 0.01 (short dashes) and a 2 — 0.04 (long dashes). 

Fig. 5: The shift in the (9-temperature of polyampholytes: ATg = Tg — Tg in units of Tg as function of the excess 
charge \z\ for different salt concentrations, eq. |6^. Only the fully charged polyampholyte case: Zj = ±1 is shown. 
The polymer concentration is 0^ = 10~ 6 A~ 3 , \v\<fr^ = 0.01 and the salt concentration: Cb = (solid curve); 
Cb = lOmM (dots); Cb = 0.1M (dashes); and Cb = 1M (long dashes). 
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